Photon Signals from Quarkyonic Matter 
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We calculate the Bremsstrahlung photon spectrum emitted from dynamically evolving quarkyonic 
matter, and compare this spectrum with that of a high chemical potential quark-gluon plasma. We 
find that while the rapidity distributions are similar, the transverse momentum distribution and the 
harmonic coefficient is markedly different in the two cases. The transverse momentum distribution 
can be fit with an exponential, but is markedly steeper than the distribution expected for the quark- 
gluon plasma, even at the lower temperatures expected in the critical point region. The V2 coefficient 
fluctuates chaotically from event to event, and even within the same event at different transverse 
momenta. The latter effect, which can be explained by the shape of quark wavefunctions within 
quarkyonic matter, might be considered as a quarkyonic matter signature. 
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The study of nuclear matter at moderate (T ~ — 
180 McV) temperature and large baryochemical potential 
(fiQ = ~ Aqcd = 250 MeV) has recently enjoyed 

new vigorous theoretical and experimental interest. 

From the experimental side, this is due to the start 
of programs specifically aimed at exploring lower energy 
collisions with the latest detector technology (THll . 

From the theoretical side, it was realized that this 
regime presents both potential for very interesting 
physics and unique challenges. The latter come from the 
breakdown of most of the techniques used to study QCD: 
lattice gauge theory presents the well-known sign prob- 
lem at finite chemical potential Several approaches 
have been invented to deal with this, but the re- 

sults are neither conclusive nor precise enough to draw 
any conclusions about the property of matter in the ex- 
perimentally interesting region. Effective field theory is 
also problematic, since the typical momentum exchange 
is ~ fj,Q, while the "fundamental scale of the theory" is 
Aqcd ~ A*q- Hence, we expect effective field theories to 
be unreliable. Indeed, while such theories generally pre- 
dict a critical point, its location on the phase diagram is 
highly sensitive to the details of the model @ . The ef- 
fect on deconfinemcnt, a non-perturbative phenomenon, 
on chiral symmetry in the critical region adds an extra 
complication which is little understood [13, [H| • 

This ambiguity leaves room for qualitatively new phe- 
nomena, and even new phases of matter, to arise. A re- 
cent proposal of this kind is quarkyonic matter fl2T4l9j . 
It is motivated by the ambiguity of defining confinement 
in a system where density ps is high enough that there is 
~ O (1) baryon per baryonic volume. The large number 
of colors (N c ) expansion [201422T ] can provide some in- 



sights on the properties of matter in this regime. While 
rough, this expansion is immune from the above issues. 

The possibility of quarkyonic matter (l2| comes from 
the asymmetry between the confinement scale in temper- 
ature and chemical potential: The amplitude of a gluon 
loop at finite temperature ~ iV 2 , while a quark- hole loop 
at finite chemical potential \i q has amplitude ~ N fN c p? q 
[12]. While at high temperature (T > A QC d ~ N°) 
confinement is broken by gluon loops alone, because 
of asymptotic freedom, at low temperature quark-hole 
loops need to overpower gluon loops. This requires 
H q ~ y/ 'N f I N c Aqcb a t one loop, and an even higher 
exponent z (1/2 < z < 1) at more than one loop [23j. As 
the number of colors might be considered "large", this 
introduces an extra scale \J N c Aqcd in momentum space 
relevant for deconfincment at finite chemical potential. 
In configuration space, however, the only relevant scale 
is the inter-particle distance, which for one baryon per 



baryonic volume [i q ~ Aqcd is always ~ N c 
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Dense matter at Aqcd < P"Q < \J N c /Nf Aqcd, with 
features of asymptotic freedom in configuration space but 
features of confinement in momentum space, at Aqcd < 
(J>q < \J N c /Nf Aqcd is known as "quarkyonic". This is 
an interesting idea, but how much of quarkyonic dynam- 
ics survives at N c = 3 and Nf = 2, 3 is an open question. 
It has long been known [U [22[ that there are significant 
qualitative differences between the N c — > 00 limit and 
N c = 3. Baryons in the N c — > 00 regime have an exci- 
tation energy <C Aqcd and strong binding (binding en- 
ergy goes as the baryon mass). The critical point for the 
nuclear liquid- gas phase transition is at T,[1q ~ Aqcd- 
Neither of these holds true in the real world. 

As argued in p3l - [25| . this indicates that the large N c 
limit is separated from the real world by a percolation- 
type phase transition. The quarkyonic matter transition 
line is therefore bound to be curved in N c as well as T, jiB 
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space, the former being accessible only on the lattice. In- 
deed, some aspects of quarkyonic matter (the "Skyrmc 
crystal" phase) are likely to be captured by the liquid- 
gas phase transition in the physical world [26|-|28( . Other 
aspects, most importantly the appearance of quark de- 
grees of freedom, could manifest themselves in our world 
[23| provided the deconfinement phase transition is far 
enough on the baryochemical potential axis. At present 
theory cannot tell us if this is the case, particularly be- 
cause of the cross-over nature of confinement (itself sepa- 
rated from the large N c limit by a critical point [29j). The 
existence of quarkyonic matter, having the properties of 
[lH is therefore a matter for experimental investigation, 
necessitating a quarkyonic matter phenomenology. 

To proceed further, we need to state explicitly what 
we mean by quarkyonic matter: We base our defini- 
tion on the reasoning in |l2j , recapping the consequences 
explored in section 4 and 5 of [2J|. Quark degrees of 
freedom make an appearance, and their interactions are 
governed by the Feynman rules of perturbative QCD. 
The equation of state at equilibrium, therefore, is simi- 
lar to that of an asymptotically free gas of quarks with 
a Fermi surface at fj, q ~ (1 — 3)Aqcd and low tempera- 
ture. Unlike "real" pQCD, however, confinement is still 
there: Baryons continue to exist, and quark wavefunc- 
tions are localized around baryons. As in the large N c 
limit, baryons are also approximately classical objects, 
well localized in position; They are also dense enough 
that there is, on average, one baryon per baryonic vol- 
ume O (Aqq D ^ or more. Hence, quark wavefunctions 

are not the asymptotically free quark wavefunctions of 
pQCD but are instead Eigenfunctions of a series of poten- 
tial wells at the location of the baryons [23| (Fig. [T]). This 
is very similar to the dynamics of a free gas of electrons 
in a metal, where atoms are classical potential wells (as 
baryons are at large N c QCD) and electrons are fermions 
weakly interacting with each-other but with wavefunc- 
tions determined by classical potentials (as quarks are 
supposed to be in quarkyonic matter). Summarizing, 
any dynamics inside quarkyonic matter will have pQCD 
interaction vertices, but incoming quark lines will pick 
up a form factor, reflecting their confinement. Unlike 
in the vacuum, this form factor will not be "localized" 
(since percolation is naturally interpreted as the der- 
ealization of quarks but will reflect the dynamics 
of all baryons of the system (23[. At a single time step 
in configuration space the quark wavefunction looks like 
(arrows indicate a 3- vector, Greek indices a 4- vector) 
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and <f>{x — Xi) arc peaks centered around the baryon lo- 
cation Xi with wavefunction width in configuration space 
~ Aq^ d , the confinement scale. We approximate (j) by 
Gaussian wavepackets 




FIG. 1. (color online) The form of the quark wavefunction in 
quarkyonic matter, in ID, as a red dashed wavy line. Baryons 
are represented by semiclassical black potentials wells. Ar- 
rows depict the motion of the baryons. 



the advantage of this ansatz is that the mean field in a 
given event can be solved analytically: The quark density 
in momentum space, assuming a baryon is a classical 
mean field of quarks, will be $ (&), where 



where <p are the baryonic quark wavefunctions 



4>(k, Xi) oc exp 
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The space coordinates at each time-step extracted 
from a UrQMD [3(1 HH simulation. Note that the config- 
uration space position of the baryon enters the wavefunc- 
tion as a phase factor, to be multiplied with momentum. 
Due to these phases, the full scattering rate will pick 
up interference terms, a crucial effect when momentum 
anisotropy is computed. We also note that UrQMD bet- 
ter captures the "classical" nature of baryons than trans- 
port models not based around molecular dynamics, where 
photon results are also available (32[. For the results 
shown here, we are using Pb-Pb collisions at y/s = 7.7 
GcV, which is well within the FAIR energy regime. We 
expect the energy and system size dependence of them 
to be weak, allowing comparisons with any 

of EHa. 

What are the observable consequences of such dynam- 
ics? Electromagnetic signals are sensitive to the earli- 
est, densest phase. Unlike, for example, hydrodynamic 
observables (that depend on the EoS), the form factors 
directly influence the final spectrum. Hence, the explo- 
ration of spectra of electromagnetic particles is an obvi- 
ous place to distinguish between quarkyonic phases and 
more conventional QGP. 

The first observable we look at is photon 
Bremsstrahhmg from quark-quark pQCD scatter- 
ing. The scattering matrix elements for this process 
were studied in [33j and are given by 



M 2 = L 2 (kx, k 2 



k 3 ,k i ,p)+L 2 (k 1 <-» k 2 ,k 3 <-> fc 4 ) (5) 



Xi) = cxp — (x 



Xi) 2 A 
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where, in terms of the fine structure constant e and the 
QCD coupling constant A. Here 

I 2 = -]e 2 X 2 N~ 2 (k 2 -k 4 )- 4 Tr [hl°hl P ] Tr [t 3 Z^Z> 
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and = 7a(fti -p)~V + ^(h+p)- 1 ^. 

We assume /iq is high enough to ignore the qg — > jq 
scattering process, to which the process considered here 
is subleading in the strong coupling constant. This as- 
sumption is consistent with the more general assump- 
tion that gluon propagators in quarkyonic matter are still 
cutoff in configuration space by confinement, and hence 
reactions involving outgoing gluons are suppressed (l2l - 
[lil . [28j . Similarly, antiquarks are neglected. 

We choose this observable because it is the simplest to 
calculate, and the one least depending on assumptions, 
such as the degree of local thcrmalization of quarkyonic 
matter, which we know little about. Local quark wave- 
functions plus well-understood pQCD matrix elements 
can be straight-forwardly integrated into a distribution 
of Bremsstrahlung photons. In contrast, for example, 
dilcpton observables (very promising as noted earlier [23| ) 
necessitate a quark and a hole distribution, where the lat- 
ter needs additional assumptions about the microscopic 
quarkyonic state to be quantitatively estimated. The 
photon rate convoluting the pQCD matrix and the wave- 
function of quarkyonic matter is then 



d7V 7 
d 3 p 



(0) 



where M. is the matrix element corresponding to the dia- 
grams Fig. [2]Eq. [5l production of a photon by the strong 
scattering of two quarks [33[ and ^ are given by Eq. [3j 
For the quarkyonic phase, Eq. [6] is calculated for each 
time step in each event, with the quark wavefunction re- 
flecting baryon location for that event. An average over 
UrQMD events is then obtained. The integral in Eq. [5] 
was done by Monte Carlo, and the results show, in partic- 
ular the fluctuation in the last figure, were verified to be 
independent of statistics, both within and across events. 

As a comparison, we also present the rate for an ex- 
panding thermalized quark-gluon plasma. The only dif- 
ference, of course, is that the incoming quark distribution 
functions are boosted-thermal, with temperature T 



*(fc)f *(k') <xS(k' - k) exp (-u^/T) . 



(7) 



Flow includes longitudinal expansion across the kine- 
matic range parametrized by longitudinal flow rapidity 
Ul , and a transverse expansion vt with anisotropy v 2 t ■ 
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( cosh(yi) 

(v + V2T cos(2</>)) cos(</>) 
(v + v 2T cos(20)) sin(</>) 
smh(y L ) 
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Confinement is incorporated by removing quarks and glu- 
ons with less momentum and virtuality than Aqcd- 

/ d 3 kl,2,3,4 -> / ^1,2,3,4*1,2,3,4 / ^1,2,3,4 



xe ((h + fc 3 ) 2 - a^ cd ) e ((fc 2 + fc 4 ) 2 - a 2 qcd ) 




FIG. 2. (color online) The pQCD process we are examining . 
The emitted photon is denoted by a wavy line, the spring-line 
is a gluon while the solid lines are quarks 



which also takes care of collinear divergences. 

In the following, we plot the experimentally interesting 
rapidity (y)and transverse momentum (j>t) distributions 
as well as harmonic distributions w.r.t. the reaction plane 
Acf> = cj) — <f>Rp. The latter is parameterized with the v n 
coefficients, of which v 2 is the most studied example 



dN- ( 
d 3 p 



diV 7 
dydptp 



1 + 2 ^2 v " (JPt, V) cos (A</>) 



(9) 



Missing in this calculation is any quark boost due to 
finite baryon momentum. Quarks are treated as essen- 
tially static. Due to the derealization of quarks, the 
consequences of finite momentum baryons are actually 
not so trivial, but they will always come with a -/V" 1 fac- 
tor, consistent with the hierarchy between "light" quarks 
and heavy baryons. In a N c — 3 world, corrections 
might be significant, although strong binding of nuclei 
in a dense medium will strengthen our argument. To 
estimate qualitatively the effect of these corrections, we 
choose a baryon at random in UrQMD, and "localize" the 
quark to that baryon, Lorentz-boosting the quark by the 
baryon's momentum. Since baryons have some flow on 
average, this boosts the flow of the quarks. This distin- 
guishes "quarkyonic" from "boosted quarkyonic" in the 
plots. While this is not a quantitative estimate, it is an 
"extreme scenario" , where a delocalized quark receives 
a "full boost" from one particular baryon. Hence, in a 
sense, it provides an upper limit to how large the bary- 
onic flow contribution can be without quarks becoming 
the actively flowing degrees of freedom. 

Similarly, this model is non-causal, since the quark 
wavefunctions adjust to baryon movement instanta- 
neously. This is another artifact of the approximation 
discussed above, fixed by ~ A^T 1 contributions, but 
not improved by the localization ansatz. For longitudi- 
nal dynamics, where typical baryon longitudinal rapidity 
yi, 1, this could be a significant issue, but for transverse 
dynamics, where baryon speeds yT <C 1, this can safely 
be ignored since quarks are much faster than baryons. 
Our most interesting results are indeed transverse. The 
results are shown in Fig. [3] for transverse momentum 
distribution, Fig. |4] for the rapidity distribution and in 
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T=150MeV,v=0.1 

-- T=150MeV,y=0.5 

■ ■ ■ Quarkyonic b-Ofm 
— ■ Quarkyonic b-8fm 

■ ■ - Boosted quarkyonic b-[)t'm 
Boosted quarkyonic b-8fm 




FIG. 3. The photon rapdity distribution for quarkyonic mat- 
ter (impact parameter b — 0, 8 fm) and thermalized QGP 
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thermalT=150MeV.v=0.1 

thermal T-150 MeV ,v=0.5 
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FIG. 4. (color online) The photon transverse momentum dis- 
tribution for quarkyonic matter and thermalized QGP 



Fig. [5] for the elliptic flow of Bremsstrahlung direct pho- 
tons, at impact parameter 6 = 0,8 fm. Normalization 
is arbitrary, as it is highly dependent on the undeter- 
mined strong coupling constant in the quark-quark scat- 
tering processes. Any determination of quarkyonic mat- 
ter would come from the shape of the distributions. 

As can be seen from Figj3l the rapidity distribution 
of the quarkyonic and thermal plasma are practically 
equivalent. In boosted-quarkyonic, a ccntrality depen- 
dent double-hump structure emerges, reflecting the non- 
negligible longitudinal flow and transparency of baryons. 
The latter can however easily be accounted in a thermal 
model with a better rapidity dependence of chemical po- 
tential. Hence, we do not expect the rapidity distribution 
to be a useful distinguisher. This is perhaps not so sur- 
prising: All variables sensitive to the difference between 
QGP and quarkyonic matter, the incoming momenta and 
virtualities, are integrated out. The resulting difference is 
well within parameter uncertainties (temperature range, 
coupling constant range, Debye screenings etc.) which 
have historically plagued this field. 

Fig. 0]shows the pt distribution for quarkyonic matter 
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FIG. 5. The photon V2 for quarkyonic matter and thermal- 
ized QGP . See text for an explanation of the main result, 
the chaotic V2 spread, which is independent of binning and 
statistics. An extrapolation from hadron V2 is also shown 



is distinctively steeper than realistic thermal curves. To 
reproduce it with a thermal curve, an unrealistically low 
mean temperature has to be used, one well below decon- 
finement, even at realistically large chemical potential. A 
steep pt distribution, therefore, can be considered as a 
possible signature for quarkyonic matter. 

We note that this steepness is natural to explain within 
quarkyonic assumptions: quarks are delocalized, and 
hence do not feel the flow of any particular hadron. 
Boosted UrQMD, unsurprisingly, is much less steep, but 
still well on the low side for temperature (T ~ 30MeV), 
less than both hadronic and partonic dynamics [3414381 ]. 

While photon elliptic flow, and spectra in general, have 
yet to be measured for the energies discussed here, pre- 
vious experience [3^, H(| suggests photon flow follows 
hadron flow closely. The latter, to a good approximation, 
is energy- independent when binned with px [H|- While 
this is something hydrodynamic and transport models 
have yet to account properly (23, |42|, we can use it to 
extrapolate (thick line in Fig. [5]). 

While, unsurprisingly, the flowing ansatz reproduces 
very well the observed trends, qualitatively and quan- 
titatively with the right choice of parameters, quarky- 
onic «2 is something qualitatively different: Overall v-i 
is compatible with zero, but with strong variations both 
event by event and within the same event in different px 
bins. This spread is approximately constant with central- 
ity. This, while completely different from anything seen 
before, is actually physically not surprising: If baryons 
are not moving, the only source of vi are the effects 
of the baryon distribution on the quark wavefunction. 

— 1/3 

The latter oscillate with a frequency ~ PtPb m an 
event- by- event dependent manner, determined by both 
the density and the quark momentum. The inter-baryon 

— 1/3 

distance p B is highly inhomogeneus, both within the 
same event and event by event. Hence, both that V2 = 
overall, and its chaotic oscillation are not surprising. 
Such behavior is not seen in cither a hadron gas or 
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an equilibrated QGP system, which is natural since V2 
in these systems is generated by flow and not by wave- 
function shape. Because hadronic transport V2 at these 
energies is so low, this does not change when the boost 
term is put in. Hence, both transverse and elliptic flow 
are expected to be not significant in quarkyonic matter 
on average, but fluctuate significantly when binned in pt- 
N^ 1 corrections corrections do not change these conclu- 
sions, so photon V2 can be developed into a signature for 
quarkyonic matter. Because the photon signal is pene- 
trating, and the UrQMD V2 signal itself is weak, we can 
be confident that the approximation we made that the 
system remains quarkyonic throughout its lifetime will 
not alter the most important signature, the photon V2 
spread in pt ■ 

In conclusion, we have calculated the rapidity, pt and 
harmonic distribution of photons from quarkyonic mat- 



ter, defined as a gas of perturbative quarks moving in 
a baryon generated classical potential. We found that 
the photon spectra is much steeper than that expected 
from a QGP or a hadron gas [HI at similar temperature, 
and x>2 oscillates in a way that has never been observed 
before, but can be understood from the shape of quark 
wavefunctions. We eagerly await experimental analyses 
capable of investigating these predictions. 
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